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. Abstract. We discuss results of our study on AM CVn binaries formed 

^vq ' with donors that never ignited He before contact. For the hrst time, we treat 

the donor's in these systems in the context of a full stellar structure evolution 

y—{ ■ theory and find that the binary's evolution can described in terms of 3 phases: 

^ ' contact, adiabatic donor expansion, and late-time donor cooling. Details of the 

ON , first and third phase are new results from this study and we focus on generally 

■ characterizing these two phases. Finally, we present our predictions for the 

' donor's light in these systems. 

o ■ 
^ : 

Introduction 

o ■ 

The AM CVn stars are a class of interacting binaries with ultracompact or- 
bital periods (P rb, that is P Q rb ^ 70 min), a white dwarf (WD) accretor, and 
a compact donor with a He dominated composition. These systems are their 
progenitor systems stellar evolutionary end-state. Thus, they provide a forensic 
record of the binary evolution processes that formed their population. They will 
also be among the gravity wave (GW) sources probed by future space borne 
interferometer missions such as LISA. Determining how observations can con- 
strain AM CVn formation processes as well as developing expectations for this 
population's distribution of GW signal frequencies and amplitudes depends on 
the details and accuracy of both our population synthesis and stellar evolution 
theory for these systems. 

The population synthesis calculations of iNelemans et aL | (|200lh illustrates 



this point. These authors considered two formations channels for AM CVn 
binaries, both involving a sequence of two common envelope events that place 
the remnant stellar cores close enough that GW emission can drive them into 
contact. The channels differ in the nature of the proto-donor prior to contact: 
either a non-degenerate He-burning star or a fully degenerate He WD. In their 
model, the formation channel determined which of two predetermined mass- 
radius (M-R) relations the donor is constrained to evolve along under mass loss. 
Donors formed from He-stars are hotter and have larger radii, i?2, at fixed mass, 
M.2- Consequently, as compared to systems with a WD donor, such systems lead 
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to larger mass transfer rates, M, at a given -P or b> evolve out to longer P or b, and 
lead to a differ ent distribution o f GW signals. 

While the iNelemans et al.l (l200lh study u t ilized a single donor M-R rela- 
tion for each formation channel. iDelove et alJ (2005) argued that even within 
each channel, there should be variations in donor properties. In particular, they 
examined the impact of variations in the donor's entropy (or, equivalently, de- 
generacy) at contact for WD channel AM CVn systems. To calculate the donor's 
evolution during the AM CVn phase, they utilized set of He-object M-R rela- 
tio ns parameterized additio nally by the donor's specific entropy, s, developed 
bv lDelove k, Bildstenl (|2003h and assumed the donor's evolve adiabatically. The 
results were to shift the distribution of WD channel AM CVn systems off the 
fully degenerate track to higher M at fixed P rb and to shift the maximum P or b 
in the population to longer periods. In particular, some of the WD channel 
systems overlapped w i th the expected locations of He-star channel systems. 

The lDelove etafl {2005) calculations have several shortcomings: the donor's 
are still constrained to evolve along predefined M-R tracks completely defined by 
the donor's s at contact. These M-R relations used were calculated assuming the 
donor's are fully convective and no treatment of the donor's thermal evolution 
was included. It is unclear whether the assumptions of adiabatic donor structure 
and evolution are valid during all portions of the AM CVn phase. These models 
also do not allow predictions of the donor's light. And, finally, wi thout full stellar 
models, accurate treatment of the important contact phase (|Nelemans et al .1 
120011 : 1 Marsh et al .Il2004h , during which time the M rate grows from to its fully 
developed rate, is not possible. 

To rectify this situation, we have calculated the evolution of WD channel 
AM CVn systems from pre-contact to late times treating, for the first time, 
the donors within a full stellar structure theory. To do so, we developed a new 
Henyey-style stellar evolution code base designed to be maximally flexible in 
terms of input physics used and in definition of the sys tem of ODEs to be solved. 
We used the pure He equation of state developed bv IWinisdoerffer fc Chabrierl 
(2005) that treats the difficult region of partial pressure ionization and covers the 
wide range of parameter space needed to follo w the donor's com plete evolution. 
We ca lculated M based on the prescriptions of lRitteii ( 19881 ) and lKolb Ritten 
(1990), allowing us to follow the binary's evolution during the phase of M turn- 
on. Below we discuss the new insights into AM CVn evolution this improved 
donor treatment provides, both in terms of the binary's evolution and the donor's 
contribution to the binary's light. 



Phases of AM CVn Evolution 

We find that AM CVn evolution can be divided into 3 phases, illustrated in 
Fig. ^ The AM CVn binary's evolution initiates when the donor is brought 
into contact and M begins increasing from zero. The first of our 3 phases, the 
M turn-on phases, consists of evolution from the start of mass transfer up to 
M maximum. During this phase, both i?2 and P or b decrease. By the time 
M reaches its maximum, Ri has begun increasing and the long-lived second 
phase where the donor evolves nearly adiabatically in response to mass loss is 
underway. Finally at late times, the donor begins cooling and contracting once 
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Figure 1 . A representative result of our AM CVn binary evolution calcula- 
tions showing the evolution of i?2 and M as a function of time. The binary's 
initial M 2 , accretor mass (Mi), central degeneracy parameter ( ip c = ep >c /kT c , 
where £p jC is the electron Fermi energy at the donor's center), and P or b at 
contact are indicated on the plot. The system's evolution can be divided into 
three phases as marked by the dotted lines: the M turn-on phase before M 
maximum, a phase of adiabatic donor evolution, and a final phase of donor 
cooling at late times. 



again, initiating the third phase of the system's evolution. The details of the first 
and third phases are new results, so we will focus here on providing a description 
of them. 

The salient features of the M turn-on phase are the R2 evolution and the 
growth rate of M. The transition from donor contraction to expansion that 
occurs in this phase can be explained in terms of the donor's initial s profile. 
As Mi decreases, the pressure, P, at fixed mass coordinate, m, in the donor 
decrease, as does density, p, there. This results in fixed mass points moving 
outward in radius under mass loss. This general trend towards expansion is 
counteracted by the fact that the mass lost also removes its own R2 contribution. 
Whether the remaining regions of the star expand enough to make up for this 
lost i?2 contribution determines whether the donor expands or contracts overall. 

At fixed m, the contribution to R2 evolution depends on (din P/dM<2) m ~ 
P c j(PM-i)-, where P c is the donor's central pressure. The R2 changes are there- 
fore dominated by the response of the outer layers. To determine if the donor 
should tend towards expansion or contraction, it is sufficient to determine (din pjdM^) 
at fixed P in the outer layers: if mass elements advected to lower P arrive with 



Deloye et al. 



a 



O 



10 



-11 



-12 



1 




_ M 2 , = 0.20, M u = 0.30 

log(*„i) = l-l. 1-5, 2.0, 2 5, 3.0, 3.5 



50 



(min) 



Figure 2. The evolution of the binary's M as a function of P or b for systems 
where the initial donor and accretor masses were 0.2M Q and 0.3ATq, respec- 
tively. The different curves correspond to different initial donor entropies, 
parameterized here in terms of ipc,i- The solid line shows an essentially fully 
degenerate object, lines above the solid line show the evolution for decreas- 
ingly degenerate donors. 



higher p than the element it replaces, R2 will contract overall. One can show 



dM J p Xp \dM- 

where T is the local temperature, V = {d\uT/d\nP) describes the background 
T profile, V' describes the actual T evolution of the advected element, and 
Xt/Xp = {dhxp/dhxT) P . 

In general, \T > 0, so that for R2 to decrease with M2, two requirements 
need to be met. The first is that xt 7^ 0; i.e., the outer layers can not be strongly 
degenerate. The second is for V — V < 0. The latter condition occurs in radia- 
tive layers under rapid mass loss (i.e. when V' ~ V a d)- The second condition 
dominates the donor's initial response; that is, initially near their surface the 
donors have a very steep radiative s profile. Once M has increased sufficiently 
so that advection is nearly adiabatically, the steep s gradient leads to increasing 
p at fixed P as much lower entropy material replaces the removed layers. This 
produces the overall R2 decrease. The contraction is eventually reversed one the 
layers with a steep s gradient are removed, exposing the transition to the much 
more shallow s profile that occurs deeper in the donor's outer layers. 

During the turn-on phase, the M growth rate depends on the evolution 
of AR = R2 — Rr,- In itially, at low M, Rl > R2 and M grows exponentially 
with AR (|Ritterlll98^ . The M maximum is set by the M needed to reverse 
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Figure 3. The evolution of the donor's L, T e g (lower set of curves) without 
irradiation and of L, T p hot (upper set of curves) for an 77 = 0.1. In all cases, 
the initial donor and accretor masses were 0.2M Q and 0.3M Q , respectively. 
The different curves correspond to the same initial donor entropies as in Fig 
|2| with higher curves having higher entropy. 



the inward evolution of orbital separation driven by GW emission. In all cases 
considered, the maximum M value is greater than the M able to be supported 
if > i?2- Therefore Rl eventually moves below the photosphere and the M 
growth rate slows since Rl '\s now probing the donor's interior p profile, not the 
profile of an isothermal atmosphere. The first change in the slope of M(t) seen 
in Figgis due to this effect. The later, rapid upturn in the M growth rate is 
related to the evolution towards R2 minimum. The increasing p of the donor's 
surface layer caused by the increasingly adiabatic advection of the underlying 
lower-s layers produces this rapid upturn in M. Less degenerate donors make 
contact at wider a, and such systems have lower M maxima since a lower M is 
needed to compensate for the weaker GW emission of these systems. 

By the time the system has passed its M maximum, the donor is expanding 
adiabatically and the binary is now in what is typically thought of as the AM 
CVn phase. In prior work on AM CVn evolution, the donor's adiabatic evolution 
has been assumed to continue indefinitely. Our more detailed calculations show 
that this is not the case and, generically, the donors eventually cool and contract. 
The reason for this is most easily described in terms of the evolution of two 
time scales: the donor's mass loss time, tm = M2/M, and the donor's thermal 
time, T t h = c' P T c M2/ L, where c' P = J cpTdm/ (TqM?) rs const, and T c is the 
donor's central temperature. Once M2 is much less than the accretor mass, 
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d In tm /d In M2 ~ 4(n — 1/3) — 1. During adiabatic expansion, dlnr t h/dlnM2 ~ 
1 + V a d, c (2 — 4n) where n = (din R,2/dn\n M2) is calculated from the donor's 
evolution. 

For typical values of n and V a d,c> ^th/^Af oc M^, where a « —5.5. During 
the early portion of the adiabatic evolution phase, tm <C T th ; but this inequal- 
ity is eventually reversed due to this large negative a value. This occurs for 
Port, f» 40 — 50 minutes. At this point, the donor begins shedding entropy 
and contracting. At late times, all evolution tracks collapse towards the fully 
degenerate track. This effect is illustrated in Fig. |2J 

The Donor's Light 

The donor's luminosity, L, and effective temperature, T e fj, evolution also corre- 
sponds with the donor's evolutionary phases. As seen in Fig. |21 initially both L 
and T e fj decrease rapidly as the system evolves toward M maximum. Once in 
the adiabatic phase, both quantities plateau at a level set by the donor's initial 
s profile. Finally, once the donor's adiabatic evolution ends at long P or b, both 
L and T e g begin decreasing more rapidly again. 

We can compare the amount of accretion light the donor sees to its own 
luminous output. The luminosity generated by the accretion flow is L acc = 
M(4>Li — <Ari)j where (f>Li, <j>R\ are the value of the binary's potential at the 
inner Lagrange point and accretors surface, respectively. At almost all times, the 
accretion light irradiating the donor dominates the donor's own luminous output. 
To treat this effect, we performed calculations where we set the donor's outer 
boundary temperature to Tp hot = T e 4 ff + r]L acc / (IQ-Kcra 2 ) where rj is the overall 
efficiency with which L acc is converted to photons that end up being reprocessed 
in the donor's atmosphere, and a is the binary orbital separation. Here we define 
T e ff in terms of the donor's intrinsic L, L = AiraR^T^. The luminosity seen from 
the donor's surface, L s is then defined by L s = ATtaR^T*^. The evolution of 
L s , Tp hot for rj = 0.1 is shown by the upper set of lines in Fig. |21 As compared 
to the non-irradiated case, L s can be increased significantly over L for a large 
portion of the AM CVn phase. Overall, with 77 = 0.1, we expect donors to have 
L s in the range of 10- 6 -10- 3 L o and T phot between 1000-5000 K for P orb < 50 
minutes. Both quantities will be higher if rj is larger. 
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